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Abstract 

We consider the exclusive B s —> <^t + process in the standard model using a constituent 
quark loop model approach together with a simple parameterization of the quark dynam- 
ics. The model allows to compute the decay form factors and therefore can give predictions 
for the decay rates, the invariant mass spectra and the asymmetries. This process is sup- 
pressed in the standard model but can be enhanced if new physics beyond the standard 
model is present, such as flavor-violating supersymmetric models. It constitutes therefore 
an interesting precision test of the standard model at forthcoming experiments. 
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I. INTRODUCTION 

We discuss the B s — > exclusive process using a Constituent Quark-Meson model (CQM) El based on quark- 

meson interactions. Quark-meson interaction vertices can be obtained by partial bosonization of a Nambu-Jona- 
Lasinio type four-quark interaction vertices for heavy and light quarks (2). In this model the transition amplitudes 
are evaluated by computing diagrams in which heavy and light mesons are attached to quark loops. Moreover, the 
light chiral symmetry restrictions and the heavy quark spin-flavour symmetry dictated by the Heavy-Quark-Effective- 
Theory (HQET) are both implemented. 

Flavor Changing Neutral Current processes (FCNC), like B s — > (/>/z + /x~ where the b quark is transformed into a 
s quark by a neutral current, are absent in the Standard Model (SM) at the tree level. This makes the effective 
strengths of such processes small. In addition, these transitions are dependent on the weak mixing angles of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix and can be suppressed also due to their proportionality to small CKM 
elements. If the quark masses in the FCNC loop diagrams are close to each other, the GIM mechanism is effective 
and this implies that FCNC charm decay transitions are very suppressed. On the other hand, FCNC beauty decays 
should be observable at LHC. These decays, known as rare B decays, are extremely interesting for the study of new 
physics. Indeed, any observed enhancement of their branching ratio could be the trace of some no-SM mechanism. 
Besides the sensitivity of rare B decays to new physics, their study is very important in the context of the SM for a 
comparative determination of the CKM matrix elements Vtb, Vtd and Vt s ||: these quantities can be measured directly 
in top quark decays. 

Here we study, in the CQM approach, the exclusive process B s — > (f>fj, + . Our aim is to provide an independent 
determination of this process in a different context with respect to that of QCD Sum Rules, where it has been studied 
first [H . Exclusive semileptonic decays are in general more complicated than inclusive modes on the theoretical ground 
since they require the determination of form factors. On the other hand they are very promising on the experimental 
side, being accessible to a large number of ongoing and future experiments. The method we apply can be used for the 
calculation of other processes like B s —> K*l + t~ or B u ^ — > K*£ + £~, Bd, s — * rj£ + £~ following the same techniques 
described in this paper. 

The CQM model has turned out to be particularly suitable for the study of heavy meson decays. Its Lagrangian 
describes the vertices (heavy meson)-(heavy quark)-(light quark) Q, and transition amplitudes are computable via 
simple constituent quark loop diagrams in which mesons appear as external legs. In the case of the process at hand, 
two interfering diagrams contribute to the transition amplitude, see Fig. 1 and Fig. 2. They correspond to the two 
physically distinct alternatives in which the FCNC effective vertex is either attached directly to the quark loop or 
via an intermediate heavy meson state (the matrix elements factorize into an hadronic and leptonic part). We will 
compute the two contributions to the process and discuss their relative weight. 



1 



II. EFFECTIVE HAMILTONIAN AND FORM FACTORS 



At the quark level, the rare semileptonic decay b — > s£ + £ can be described in terms of the effective Hamiltonian 
obtained by integrating out the top quark and W ± bosons: 



H cS = -4 



„ 10 
^ 2 *=i 



(1) 



We will use the Wilson coefficients calculated in the naive dimensional regularization scheme || (see Table I). The 
numerical values used in the calculations are given in Table II. The Hamiltonian (Q) leads to the following free quark 
decay amplitude [Q: 



M(b^s£ + £-) 



Gfol 
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(2) 



where Ci(p = m&) are the Wilson coefficients which act as effective coupling constants in the 4- Fermi formulation 
of the interactions. Here C| ff = CV — C5/3 — Cg. Short-distance Wilson coefficients are redefined in such a way to 
incorporate certain long-distance effects from the matrix elements of four-quark operators Oi with 1 < i < 6. Cg ff , 
the Wilson coefficient of Og — e 2 /167r 2 (s7 M L6)(i?7 A1 ^), is defined in terms of these matrix elements in the Appendix. In 
order to compute the (<fi£ + £~\A4(b — > s£ + £~)\B S ) 1 we need the following form factors parameterization for the V — A 



terms: 



< He,p)\s^(l - T o)b\B s {p B ) > = W{ f ] e^ af} e* v p%p 

m B + m<i> 

- i^{m B + m^Axiq 2 ) 

iPB+p) »A 2 (e) 



+ i(e* • q) 
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where q = p B — p and: 



with the condition: 



m B + 
2m d> 



A^q 2 



m B - m 4 
2ms 



A 2 (q 2 



MO) = A>(o). 

The form factors needed for the magnetic term in (^) are defined by: 

< ^{e,p)\sa^q u {\ + j 5 )b\B s (p B ) > = ie lxva0 e u p' B p^2T 1 (q 2 ) 

+ T 2 (q 2 ) [e M (m| - m\) - (e ■ p B )(p B 

T 3 (q 2 ) (e-p B )q^ 



(m| 



(PB + p) A 



with: 



Ti(0) =T 2 (0). 



(4) 



(5) 



(6) 



(7) 



These seven form factors can be calculated with the aid of a Constituent-Quark-Meson (CQM) model in which the 
interactions B s — <j>— (V, A, T)current are modeled by effective constituent quark loop diagrams. The <j> is considered 
to have a photon-like interaction with the light s degrees of freedom described by: 



C 



ig^ssSjfj.sq 



(8) 



2 



where g^, ss = the factor of 3 coming from the charge of the s quark, g^ is extracted from the measured 

electronic width of the <j> via T((f> — > e + e _ ) = (47ra 2 )/3(l/<7?)m0. This way of modeling the interaction derives from 
the Vector Meson Dominance (VMD). Vector meson dominance can be obtained from the effective four-quark theory 
of the Nambu-Jona-Lasinio type for light quarks once electro-weak interactions are added. One can show f| that 
the coupling of the electromagnetic field with quarks is turned into a direct coupling of photons and neutral 
vector mesons in the effective theory, and this reproduces the VMD term. Eq. (||) can therefore be interpreted as a 
quark-quark-light-meson vertex of this kind. 

As stated above, there are two kind of contributions to the form factors, depicted respectively in Fig. 1 and Fig. 
2. In the first one the current is directly attached to the loop of quarks. In the second, there is a intermediate state 
between the current and the B s cj> system. This intermediate state is a l - or 1 + heavy meson with a s valence quark. 
The Feynman rules for computing these diagrams have been discussed in [jlj and the extension of the model to the 
strange quark sector has been developed in H]. Consider for example the direct diagram in Fig. 1. The model allows 
to extract the direct contributions to the form factors V D , Af, .., Tf, .. through the calculation of the loop integral: 

/7 i^L r S .^ Tr [ (7 • 1 + m ^ssl ■e){ 1 -{l + p)+ m)(V, A, T)^ 75 ] 

V ffmB T6W {P ~ m?)((l + p) 2 - m 2 ){v ■ I + A H ) ' ( ' 

where V,A,T denote respectively Vector, Axial- Vector and Tensor (T ~ cr M „(l + 75)) currents, I is the momentum 
running in the loop, p — m^v' is the momentum of the <f>, to the constituent quark mass of the strange quark (we have 
m = 510 MeV), v is the four velocity of the incoming heavy (CP) meson; the heavy quark propagator and the heavy 
meson field expressions from Heavy Quark Effective Theory (HQET) are taken into account. The constant Ah is 
defined as the m^ - m& difference (between the mass of the heavy meson and the constituent heavy quark mass) and 
represents, together with the cutoffs, the basic input parameter of the model. Following [Q, here we assume Ah = 0.6 
GeV. Zh is the coupling constant of the heavy meson field H (using the notations of HQET where H represents the 
(Q _ ,l _ ) heavy meson multiplet) with the quark propagators. Integrals like (|^) are obviously divergent. We define 
the model with the proper time regularization procedure which consists in exponentiating the light quark propagators 
in the Euclidean space, and introducing ultraviolet (A) and infrared (/i) cutoffs in the proper time variable. In our 
numerical analysis A = 1.25 GeV and /1 — 0.51 GeV. The results are quite stable against 10 — 15% variations of 
the cutoff values. We again refer to Jl],f7| for a discussion on the determination and the physical meaning of these 
parameters. 



A. Direct form factors 



The CQM expressions for the contributions to the form factors, derived by the direct diagram calculations with V 
and A currents (see Fig. 1), are the following: 
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(10) 
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where: 
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and: 



" = — V — > (14) 

9 9 9 

m n — q ~ mi 

ri = - J M| (15) 

The definitions of the functions Z. fl; are the following: 

z _ iNo r cg dH 

16~* J (£ 2 - m 2 )((l + p) 2 - m 2 )(v ■ £ + Ax + ie) 1 ' 

Z^ = Cl^ + ft 2 ^ (17) 

Z^ = n 3 g^ + tt^v" + {l 5 v ,fl v' v + n 6 [vV + v'»v v ] , (18) 

where Z^ and Z^ v mean that in the Z integral we are considering £^ and t^l" numerators, v is the four velocity of 
the incoming meson, v' the four velocity of the (f>. In the computations, v and v' are related by v ■ v' — (Ai — A 2 )/m0 
where Ai is v ■ k, k being the residual momentum of the heavy quark. The explicit expressions for the integrals Z , 
are algebraic combinations of some Ii integrals which are given in the Appendix. They are in general functions of the 
two parameters, Ai and A 2 . In the case of the direct diagrams in Fig. 1, A 2 = Ai — m^w, where uj has been written 
in 



(14) 



B. Polar form factors 



The polar contributions to the form factors come from those CQM diagrams in which the weak current is coupled 
to B s through an heavy meson intermediate state, see Fig. 2. The form factor will then have a typical polar behavior: 



F p (0) 
1 — g 2 /mp ' 



^(<r) = , tt^i (1 9 ) 



where mp is the mass of the intermediate virtual heavy meson state. Pole masses are given in Table IV. This behavior 
is certainly valid near the pole; we will assume that it is valid all over the q 2 range that we want to explore, i.e., also 
for small q 2 values. We find: 

V*(0) = -V2 9V XF^±^ (20) 

m B B 

Af(0)= ^ 9V p -(C ~2p.ru,) (21) 

Af (0) = -^ 9v »F ^ {m 2 B ° +m *\ (22) 

where u = m B J(2m^), while g v = m^/f^ M. By B** we mean the 1 + state of the 5 multiplct (0 + , 1+) of HQET. 
The mass of B** is taken by || to be to^** = 5.76 GeV. At present there are no experimental information about this 
state. 

As for A^(i7 2 ), we have to impose the condition (||); a choice is: 

A o(l 2 )=^(0) + 9 vf3F J / (23) 

771.0 yz?7773 3 m B 3 ~~ " 

F and F + are the lepton decay constant of the H and S HQET multiplets || : 

(VAC 1 97^75 Q | H) = i^/rTT^r^F (24) 
{YAC\q^Q\S) =iy/ra^r"F + , (25) 

where we use r = v (the heavy quark velocity) for H = 0~ , S = + and r — e (the polarization vector of the heavy 
meson) for H = X~ , S = 1 + . One finds: 
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F = 2^Z H (h + (A H + m)I 3 (A H )) 
F+ = l4z~ s (h + (A s - m)I 3 (A s )). 

The numerical table connecting Ah and As values has been discussed in 0. We notice that: 



(26) 
(27) 



fB. 



and numerically, we find: 



f B . = 180±il MeV, 



(28) 



(29) 



which is in reasonable agreement with the result from QCD sum rules g and from lattice [ [To| , according to which 
fg, — 190 MeV. The theoretical error in the determination of f Bs comes from varying the Ai = Ah parameter in 
the range of values Ah = 0.5,0.6,0.7 GeV. 

The CQM explicit expressions for the strong constants A, (3, /i, £, parameterizing the strong couplings HHcj) and 
HScj) according to the interaction Lagrangians discussed in S, are given by: 



X=^P^Z H (-n 1 +mZ) 
V2gv 

= y/2ii2Lz H [2mQ.i + mAl 2 + 2Q 3 + n A + fl 5 ~ m 2 Z\. 
9v 



(30) 
(31) 



Here the functions Z, fL are computed with Aj = A 2 = Ah, x = m$, u> = m^j '(2m Bg ) where one takes the first 
1/toq correction to lo — 0. Moreover: 
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\fZ H Z s (m^ 2 + 2ft 3 + ^ 4 + Q 5 - m 2 Z) 



(32) 
(33) 



where Ai = Ah, A 2 = As, x = and lo — (Ai — A 2 )/to^. The suffix S indicates the (0 + , 1 + ) multiplet of HQET. 



C. Direct tensor form factors 



Let us now turn to the Ti form factors. The contributions to Ti coming from the direct diagrams in Fig. 1 are 
labeled by TP. A calculation of the loop integral (||), in which we retain the T current of the electromagnetic penguin 
operator, allows to extract the Tf by comparison with (^|). We obtain the following results: 
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(35) 



(36) 



where: 
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J = (37) 

2 9 9 

^ = _^ — * (38) 

P= Bs / - , (39) 
and the following consistency condition is satisfied: 

T/V = 0) = T 2 D (q 2 = 0). (40) 

D. Polar tensor form factors 

The calculation of polar contributions to the Tj form factors follows the computation in || . In the latter reference 
a different parameterization of the tensor current matrix element is used (and also a different definition of cr M „ from 
the one adopted in this paper; namely the <7 M „ is defined without an overall i): 

(<j>(p, e)|s<7 M „(l + ~f 5 )b\B s (p B )} = i[{pB^v - PBu^^ ~ ie l _ ll ,\ a p B e' y )A(q 2 ) 

+ ip^u - Pvt/i - ie tJ .uXcrP X e' 7 )B(q 2 ) 

+ (e • Pb)(3PbhPv - PBvPti - ie^XaPBP^Hil 2 )]- ( 41 ) 

The relations between the form factors A(q 2 ), B(q 2 ), H(q 2 ) and our form factors Ti(q 2 ), T2(q 2 ) and T^{q 2 ) are given 
by: 

T 1 (q 2 ) = - t -[A(q 2 ) + B(q 2 )] (42) 
Uq 2 ) = ~\ [A(q 2 ) + B(q 2 )] - j [A{q 2 ) B{q 2 )} -^^ (43) 

T 3 (q 2 ) = j [A(g 2 ) - B(g 2 )] + ^(<Z 2 ) (m| s - m 2 ) . (44) 

Again, the condition T^O) = T 2 (0) is manifestly satisfied. Now, the polar contributions from the 1~ and 1 + interme- 
diate states have been computed in || with the following results; if the 1~ state is taken into account: 

_ (45) 
(mj, - q z )^/m Bs 

-a^FXqvm' 2 

B <" " sr? (46) 

, . ; 7 2v y A ^ . (47) 



{jTip - q 2 )y/m B 



If instead we consider the 1 + contribution 



■ 4 <' » = (48) 

5(g 2 ) = (49) 



2 -i2y/lmB~ s F+fig v 

H{q ) = — — ^ . (50) 

[Trip — q A )mB. 



where mp is the mass of the intermediate polar state. Let us consider the contribution due to the 1 state using the 
results for A(q 2 ) and B(q 2 ) obtained in ||. We find: 

Tf(, 2 =0) = --^=T/V = 0), (51) 
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neglecting a term — — -. The contribution to these form factors due to the 1 + state is sub-leading, being Ofl/ms). 

The form factor T 3 P (0) has the same structure ( |5l| ) of Tf 2 (0) for the 1 contribution. The sub-leading contribution 
from 1 + is instead: 

If (0) = (^f(P + MC ~ 2^))-L + ^^^(^)) • (52) 

We do not include the sub-leading contributions in the numerical analysis since they turn out to be very small 
corrections, certainly below the theoretical error induced by the model itself (varying e.g. the parameter Ah in the 
range 0.5,0.6,0.7 GeV). 



E. Numerical results for the form factors 



The form factors used in the branching ratio calculation are obtained adding up polar and direct contributions: 



W)=^(g a )+^(9 a ) 



(53) 



The q 2 dependence of the form factors is given in Fig. 3 and 4. The dependence on the model parameter Ah is less 
than 10%. Fig. 3 is similar to what obtained for the decay B — > plv (note however that Fig. 3 of jll| contains a 
misprint and A\ and A 2 are interchanged) . As in the B — > p decay the V form factor results in the model from a large 
cancellation between the direct and polar term. Therefore the result for the V form factor has a large incertitude. 

In order to compare with results from other approaches we consider the following two parameterizations of the form 
factors: 



F(q 2 



ap 



(*) + »' G 



(54) 



and 



F{q 2 ) = F(0)exp 



ci- 



C-2 



C 3 



' B, 



(55) 



The coefficients ap, bp, c\, C2 and C3 are given in Table III. Results can be compared with those obtained from QCD 
sum-rules (see Table III and Fig. 4 of that paper). In order to allow an easier comparison we write their F(0)sr in 
Table III. As already stated our V form factor is affected by a large incertitude. Concerning Aq, A\, A2, their value 
at q 2 — is different in the two models but their shape as a function of q 2 is quite similar. The tensor form factors 
have a more pronounced pole-like behaviour in the QCD sum-rule calculation than in the present one. 

III. RELATIONS BETWEEN THE FORM FACTORS 

Relations between the form factors can be established using the equations of motion for the heavy quarks or taking 
limits of the general expressions calculated in the preceding sections. They are useful to link different decay processes 
and as a cross-check of the calculations. 



A. Semileptonic and tensor currents 



The equations of motion of the heavy quark implies 

1+/ 



b = b 



(56) 



and in the b rest frame one has 



7 °6 = b 



(57) 
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which can be used to relate vector and tensor currents [Q 



(58) 



Therefore the form factors 7\, T2, T3 of eq. (|6|) can be related to those describing the weak semileptonic transition 
B — > of eq. (||) or £? — > p, using SU (3) symmetry. 

Using (m), the form factors Ti, T2, T3 are related to the form factors V, A\ and A2 as follows: 



ri(g 2 ) 

T2(q 2 ) 



1 



2m b b 



m B 3 ~ m l + Q 2 



m,B s + m<t, 

1 r n? 2 ) 



V(g 2 ) - (m Ss + m^)A 1 (q 2 ) 



;{; 



2m Bg (m| s - m 2 ) I mB, + to^ 
Ai(q 2 )(m Bs + m <il )(m Bs ~ m^ + <? 2 ) j 

- 1 —— Av{q 2 ){m 2 Bs + 3m 2 - q 2 ) + A 1 {q 2 ){m Bs + m^f 



q 2 



■(m| s - m|)(m| < - to 2 . + g 2 )|. 



(59) 



(60) 



(61) 



In a similar way the other equivalent parameterization in terms of A, B and H is related to the form factors V, A\ 
and A 2 : 



A(q 2 ) = 1 
B(q 2 ) = l 



.. f q 2 - M 2 - m K , 



V(q 2 



M B + m K * 



Mb + rn K * 
Mb 



Al(q 2 ) 



2M f 



Mb + mK- 



■V{q 2 



2 1 



H{q Z ) = TT- 



V(q 2 ) 



1 q 2 + M B - m 2 K , 



M B {M B +m K , 2q 2 M B +m K * 



A 2 (q 2 



(62) 
(63) 
(64) 



These relations are strictly valid only for q 2 



B. Final Hadron Large Energy Limit 

We examine a particular limit for the B — > <fi semileptonic form factors, the one of heavy mass for the initial 
meson and of large energy for the final one (LEET) . The expressions of the form factors simplify in the limit and for 
B — > Vlv, they reduce only to two independent functions Jl3| ]. The four-momentum of the heavy meson is written as 
p = Mhv in terms of the mass and the velocity of the heavy meson. The four-momentum of the light vector meson is 
written as p' = En where E — v ■ p' is the energy of the light meson and n is a four- vector defined by v ■ n = 1, n 2 = 0. 
The relation between q 2 and E is: 



The large energy limit is defined as 



q z = Mjj - 2M H E + m 2 v (65) 



A Q CD,m v « M H ,E (66) 



keeping v and n fixed and my is, in our case, the mass of the <fi. The relations between the form factors appearing in 
the LEET limit constitute a powerful theoretical cross-check of the formulas derived in the previous sections. Note 
that in our model the polar diagram of Fig. 2 is sub-leading with respect to the direct diagram of Fig. 1 in the LEET 
limit; therefore only the direct part of the form factors contribute to the expression of the LEET form factors. In 
agreement with the results obtained in |14| we find the following result: 



Mq ) = { 1 ~ M^E ) C ll( Mff '^ + J^ C±{Mh ' E) (67) 

M 2 f C±(M H ,E) (68) 
M H + m v 
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my \ 

my 



a(M H ,E)-^C\\(M H ,E) 



The explicit expressions for £| | and £j_ are as follows Q : 



(69) 
(70) 



C\\(M H ,E) = 



\/M H Z H m\ 



2Ef v 
+ 4:A H m v Z 



' V 2 
E 



(71) 



U(M H ,E) = 



y/M H Z H mj 
2Ef v 



[h(A H )+m 2 v Z] 



E 



(72) 



where terms proportional to the constituent light quark mass m have been neglected. Note that to obtain the correct 
results for the <fi meson one has to replace: 



fv 



V6 



3cos(39.4°) 



(73) 



in order to take into account the constituent quark structure of the <fi meson |l5| ]. The angle 39.4° in (|7^) is the u>—(f> 
mixing. It is interesting to note that in LEET one can also relate the form factor 2\, T 2 ancl ^3 to the semileptonic 
ones and to the £j_ and £j| form factors of the LEET limit p3| : 



T 1 (q 2 ) = C±(M H ,E), 

„2 



1 



M 2 H -m^ 



(±(M H ,E) 



(74) 
(75) 

(76) 



£j_ and C|| obtained in this way agree with those of ( |7l| . |72| 



The dilepton invariant mass spectrum for the decay B 



IV. DECAY DISTRIBUTION AND ASYMMETRY 

&£ + £~ can be written in terms of the adimensional 



masses q 2 = q 2 /m 2 B , me = mijmB s and m^ = m<f,/mB s [fL6| 



dr _ G% a 2 m% e 
dq 2 



2 10 7T 5 



2-2 ^ 



1 



4m| 
1 

2m 2 



|cr(»? - j ) + Isr h - y + 4 ™'( 2 + 2 ™0 - r) 



Re(bc*)( V - !-)(! - m 2 - f ) + Re(/ ff *)(fa - s -) (I - m 2 - f) + Amfo) 



2^r) [Mfh*) ~ Re(gh*)(l - m 2 )] + ^q 2 V \h\ 2 



(77) 



where 
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and 




(79) 

The functions a to h contain the form factors dependence and the Wilson coefficients (see the Appendix). The 
invariant muon mass distribution for the decay B s — > <j)/j, + fj,~ is given in Fig. 5. Integrating the differential decay 
distribution (|77| ) allows to compute the branching fraction for the B s — > (f>/i + fi~ decay: 

B(B S -► <M + aO = 8.8 x 10~ 5 . (80) 

Note however that this number is model dependent not only due to the form factors but also to the way cc resonances 



are taken into account (see formula (124) in the Appendix). For a comparison we calculated the same branching 



fraction excluding the effect of the cc resonances: 

B(B S -> ^/Onon-rcsonant - 2.5 X lO" 6 . (81) 



This amounts to use Eq. (121) instead of Eq. (124) for the calculation. Finally in order to have a more realistic 



estimate of the branching ratio we use the complete expression (124) but exclude the cc resonance regions 2.9-3.3 
GeV and 3.6-3.8 GeV from the integration as in |l7j| : 

B{B S -> <M + AOcx P -iikc = 1-9 x 10~ 6 . (82) 
The differential forward-backward asymmetry is given by |l8] ] 

dA FB r«n d< _ d 2 r + r° d< _ d 2 r ^ 

dq 2 J Q d^dq 2 J d<,dq 2 

For B s — > <f)£ + £~ decays we obtain: 

^jSr = ^ \^V tb \ 2 q\(ff [Re(6e*) +Re(af*)] . (84) 

In Fig. 6 we plot the differential forward-backward asymmetry normalized to the differential decay rate: 

dApB dAFB i dT 



l^i ■ (85) 



The position of the zero q 2 is given by 



dq dq 2 dq 



Note that in LEET the ratios T\jV and T^JAy are simple functions of rh^ and q 2 with no hadronic uncertainties, 
as can be seen from formulas (|67|-[76|). For a detailed study including radiative corrections see p3|. Decays such as 
B s — * (f)£ + £~ involve the Wilson coefficients Cf , C| ff and Cio. In extensions of the SM they can assume rather 
different values from those expected in SM. In particular the position of the zero in the forward-backward asymmetry 
is a measure of C7 /Re(Cg) which, as shown above, depends on form factor ratios. This decreases the model dependence 
of this number. Moreover the sign of C7 can be opposite to the SM one in beyond-SM scenarios. All these elements 
explain the relevance of the experimental study of the forward-backward asymmetry and the need of form factors 
computations, despite of their model-dependent nature. 
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V. CONCLUSIONS 



The exclusive process B s — > </>/i + /i~ belongs to a set of processes, like B — ► K* , B — > A"/j + /j~, B s — > r\y + y^T 1 
which will be accurately studied at B-factories. In this paper we have examined _B S — > in the framework of 

a Constituent Quark Meson Model. The meson is coupled using the VMD hypothesis. The model has extensively 
been tested in a number of exclusive processes |l]] . It provides results in good agreement with experimental data and 
with those obtained using QCD Sum Rules. In order to have a better understanding and a check of the form factors, 
we have considered the LEET limit of the B s — > (f>[i + fi~ decay obtaining consistency among V,A and T direct form 
factors. We have studied the decay distribution and the forward-backward asymmetry. As shown in |nj , CQM offers 
a versatile calculation framework also to study more exotic processes involving higher spin meson states. This aspect 
of the model will be very useful for the study of higher spin B and B s rare meson decays as soon as data on these 
states will become available. 

APPENDIX 
A. Integrals 

We list the explicit expressions for the integrals used in the text: 

o _ H-x/2) - I 3 (x/2) + (4f3(Ai) - J 3 (A 2 )] [Ax-wx/2\Z 

ni 2x(l-w*) 1-c 2 (87) 



o _ -J 3 (Ai) + J 3 (A 2 ) - uM-x/2) - h(x/2)} [x/2-A^Z 

n ' 2 - 2x(i-u,*) i - ^ 



K x 2uK i - K 2 -K 3 

"■ = T + 2(1-^) ^ 
o - Rl j_ 3K 2 -6loK 4 + K 3 (2lo 2 + 1) 

^ = 2(1-0;*) + 2(1 -c 2 ) 2 (90) 

-K x 3K 3 - 6ljK 4 + K 2 {2uj 2 + 1) 

" 5 " 2(1-^) + 2(1 -c 2 ) 2 (91) 

A> 2K 4 (2co 2 + 1) - 3co(^ 2 + K 3 ) {Q0 . 

n " - 2(1 ~u 2 ) + 2(1-^)2 ^ 

z 



I6n 4 J (k 2 - m 2 )[(k + q) 2 - m 2 ](v ■ k + A x + ie) 
I 5 (A 1 ,x/2,u)~ I 5 (A 2 ,-x/2,lo) 



(93) 



2x 

where x = m^, to, A±, A 2 , are specified in the text and K^s are given by: 

K x =m 2 Z-I 3 {A 2 ) (94) 

K 2 = A\Z- ^/2)- h(- X /2) x + 
Ax 

V ^ 7 i J 3( A l) ~ 3/ 3(A 2 ) U) . , 

A 3 = — ZH h — [Aii 3 (Ai) - A 2 i 3 (A 2 )J (96) 

xA lv A 1 [/ 3 (A 1 )^J 3 (A 2 )] , I 3 (x/2) - I 3 (-x/2) 
K ± = — Z+ £ + 4 ' (9?) 

expressed in terms of the Ii integrals, regularized using the Schwinger's proper time regularization method: 

1 I6ir 4 J (k 2 - m 2 ) 16^ 2 \ ' A 2 ' y? J 1 j 

7 3 (A): ' ' 



16tt 4 7 (fc 2 - m 2 )(v ■ k + A + ie) 
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N c f 1 ^ 2 ds 



/ B (Ai,A 2 ,w) = 



16tt 3 / 2 7 1/A 2 a 3 / 2 



e-^ 2 -^ (l + erf(Ax/i)) 
d 4 fc 



16tt 4 7 (fc 2 - m 2 )(w • fc + Ai + ie)(v' ■ k + A 2 + ie) 
1 



l! dx i 



16tt 2 



+ 2x 2 (l-w) + 2x(w-l) 

- / a e - s (™ 2 - 2 ) S -V 2 (i + erf^v^)) + 

hi 

ds e - s(m2 - 2 ^ a' 1 



16tt 3 /- j 1/A 2 

6 r 1 ^ 



/A 2 



Here we have used the definitions: 



r(a, xo,£i) = / dt e~* t 

«0 



t j.Q-1 



erf(z) 



ct(x, Ai, A 2 ,w) 



Jo 



dxe x 

Ai (1 -x) + A 2 x 



VT+2 (w-1) x + 2 (1-uj) x 2 



B. Integrals in the LEET limit 

We list the expressions for the integrals used to compute the form factors in the LEET limit: 

^ EET = Mx/2) h(-x/2) + 2xA,Z] 

^3 EET = -g^ Hh(x/2) h(-x/2)) + 2A 1 (7 3 (A 1 ) + 2x 2 Z)] 

^ EET = ^[/ 3 (x/2)-/ 3 (-x/2)] 

"e EET = ^ [x(/ 3 (x/2) - h{-x/2)) + 2A 1 (7 3 (A 1 ) + 4x 2 Z)] , 
where Ai = Ah, x — m<f,. 



C. Auxiliary functions 

The auxiliary functions introduced in the formula for the invariant mass spectrum for the B s 
defined as 



f) = T^—cf(f)v(f) + ^±^cf Tl (f) 



b(q 2 ) = (l + m 4 



1 - raj 



^2 

(1 - m )C 9 off (<f )A 2 (<f ) + 2(m b - m.)Cf f T 3 («f ) + l—^T 2 (q 2 ) 
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d(f) 


= ^[Cf(f)((l + m <t> )A 1 (q 2 ) 


- (1 - m )A 2 ( ( z 2 ) - 2m^ (<? 2 )) - 2(m h - m s )CfT 3 (q 2 )] , 


(114) 




= -r^—C 10 V(q 2 ) , 

1 + nicf, 




(115) 


m 2 ) 


= (l + m^CwAxif) , 




(116) 


9(f) 


= —^-C W A 2 {?) , 

1 + m</, 




(117) 


h(f) 




-TO )A 2 (9 2 )-2m A o (g 2 )] . 


(118) 



The Wilson-coefficients Ci are those of || (see Table I). The effective coefficients Cf s are defined as follows 

Ccff / I S~1 /Q f * 

Cf (<f) = C 9 + T(<f ) , 



(119) 
(120) 



where T(g 2 ) is matrix elements of the four-quark operators (for a detailed discussion on the perturbative and non- 
perturbative contributions see Jl6| ) . A perturbative calculation yields [^(J : 

TpertW 2 ) = H(z, f) (3d + C 2 + 3C 3 + C 4 + 3C 5 + C 6 ) 



N(l, q 2 ) (4C 3 + 4C 4 + 3C 5 + C 6 ) - -N(0, <f) (C 3 + 3C 4 ) 



where z = m c /mb and 
N(z,g a ) = 



8 , m b 
— In 

9 n 



+ -(3C 3 + C 4 + 3C 5 + C 6 ) , 



-InzH h -x - -(2 + x) 1 - x 17 

9 27 9 9 V Jl 1 



In 



l-a-l 

2 arctan - 



x < 1 
x > 1 



H(0,, 2 ) = --ln- + ---ln ? 2 + r7 r 



(121) 

(122) 
(123) 



and x = Az 2 /q 2 . In order to incorporate the non-perturbative contributions to T we follow the phenomenological 
prescription of fl8f| where cc resonance contributions from J/ip,%p' 7 . . . parametrized in the form of a Breit-Wigner, 
are added to the perturbative result: 



T(q 2 ) 



pert 



(q 2 ) + ^ {3d + C 2 + 3C 3 + C 4 + 3C* 5 + C 8 ) 



E < 

Vi=tl!(ls),...,ip(6s) 



r(Vi -^£ + £-)m Vi 
ray 2 — q 2 — imy i Y\ 



(124) 



The numerical values used for the masses, widths and branching fractions of the 1 charmonium resonances are 
given in Table V (data from [PH). The factors K% correct for the naive factorization approximation. K\ is calculated 
comparing the experimental B s — > J/tp(lS)<f) branching fraction to the one predicted by our calculation using the 
J/ip — ► fi + fj," experimental branching fraction: 



yexp • 



(125) 



The q 2 integration range for the calculated branching is (mj/^ — Tj/^) 2 < q 2 < (mj/^, + Tj/^) 2 , which is a 2Tj/^ 
interval around the J/ip resonance. The result is that a K\ = 1.36 is needed to correct for the factorization result. By 
taking 4r »j for the integration interval around the J/ip resonance only changes the branching ratio B(B S — > J /ipcp) 
from 9 to 11 x 10~ 4 which is within the experimental error. For the other Ki values we take again 1.36 as no 
experimental values are known for the higher 1 charmonium resonances. Note that for the ST/(3) related decay 
B — > J/ipK* the k factor is estimated to be ~ 2.3 ||| using inclusive B — ► X s l + l~ data. However using exclusive 
data a smaller k ~ 1.6 is obtained fllql. 
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TABLES 





c 2 


c 3 


c 4 


c 5 


C 6 


Cf 


0> 


Cio 


-0.248 


+1.107 


+0.011 


-0.026 


+0.007 


-0.031 


-0.313 


+4.344 


-4.669 



TABLE I. Wilson coefficients used in the numerical calculations. Cf = CV — Cs/3 — Co. 
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mw 


on A 1 ^\ T 

oU.4i bev 


mz 


yi. lou / uc v 


sin 6\y 


n 99^^ 


/i yiix cutoii ) 


u.oi Kjtev 


A ( U V CUtOII ) 


1 9£ T 

i.zo Lxev 


A 


u.o Lxev 


m s (constituent) 


U.ol Uev 




1 09 PnV 

± .UZ ^orc V 


u 


u.z4yi Liev 


m c 


1 9^ r^A^ 


nib 


4.0 beV 


mt 


1 l u.o vjc V 


fj, (scale) 


m 6 


1/aw 


129 


a s (m z ) 


0.119 


\v t *v tb \ 


0.04022 


\v t *v tb \/\v cb \ 


1 



TABLE II. Values of the parameters used 
in the numerical calculations. 





r 




A 2 


A 


Xi 


T 2 


T 3 


F(0) 


0.20 


0.59 


0.73 


0.29 


0.42 


0.42 


0.36 


F(0) S R 


0.43 


0.30 


0.26 


0.38 


0.35 


0.35 


0.25 


Op 


+0.65 


-0.11 


+0.78 


+2.70 


+0.78 


+0.85 


+0.62 


b F 


+0.96 


+0.49 


-0.52 


+2.30 


+0.07 


+12.9 


-0.88 


Cl 


+0.41 


-0.19 


+ 1.03 


+3.65 


+0.75 


-1.89 


+0.48 


C2 


+0.73 


+0.01 


-0.85 


-2.15 


+0.34 


+8.41 


+1.43 


C3 


-2.47 


-0.61 


+3.47 


+1.60 





-25.3 


+1.22 



TABLE III. Values of the parameters of Eq. ( |5l| ) and Eq. ( |55| ) for the B — ► <f> form factors 
in the constituent quark model using the central value A^j = 0.6 GeV. Varying A^j by ±100 
MeV gives a variation of 10% or less in the values of the form factors. i ? (0)sR is the form factor 
value in zero from the QCD sum-rules. 
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Form Factor 



Pole mass [GeV] 



V 5.416 

A 1 5.75 

A 2 5.75 

A a 5.75 

Ti 5.416 

T 2 5.416 

T 3 5.416 



TABLE 


IV. Values of the 


mass 




poles for the polar form factors. 








M MnS) [GeV] 


I>(„s) [GeV] 


Br(^(n5) It) 


j/1> 


3.097 


8.7 x 10~ 5 


5.88 x 10~ 2 


i,{2S) 


3.686 


2.77 x 10~ 4 


1.03 x 10~ 2 




3.77 


2.36 x 10~ 2 


1.1 x 10~ 5 


i,{AS) 


4.04 


5.2 x 10~ 2 


1.4 x 10~ 5 


i>(5S) 


4.16 


7.8 x 10~ 2 


1.0 x 10~ 5 


^(65) 


4.42 


4.3 x 10~ 2 


1.1 x 10~ 5 



TABLE V. Masses, widths and leptonic branching ratios of the 1~~ cc 
resonances. For J/tp and ip{2S) the branching ratio is the one to For 
the other resonances is ip(nS) — > e + e~. 
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FIGURES 



1.75 
1.5 
1.25 




d> q2 Gev2 

Y FIG. 5. Dimuon invariant 

FIG. 1. Direct diagram with V,A,T currents changing b -> s. mass distribution of the B s -> tj>n+ fjT decay. The vertical axis 
The momentum carried by the heavy quark is ttiqv + i + k where j s j n nn n s Q f 10 — 5 GeV -2 
k, the residual momentum k ~ Aqcjj, is due to the interaction of 
the heavy quark with the light degrees of freedom. 



l 

0.75 




T 2 4 6 8 10 12 14 

q 2 GeV 2 

FIG. 2. Polar diagram from 1" and 1+ intermediate states. FIG. 6. Normalized forward-backward asymmetry of the 

B s — > (j>fi+fi~ decay. The two peaks at high q 2 values are due 
to the J/ip and ip(2S) resonances. 

3 6 9 12 15 




' 1.5 



q 2 GeV 2 

FIG. 3. q 2 dependence of the semileptonic B s — > <j> form fac- 
tors V (dash-dotted), Ai (large dashes), A2 (small dashes) and Aq 
(continuous line). 
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FIG. 4. q 2 dependence of the tensor B s — » <f> form factors Ti 
(continuous line), T2 (large dashes) and T3 (small dashes). 
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